Physical Review Letters, In Press, arXiv:physics/0406063 



o 
o 

(N 

Or 



i 



c*2 

11; 

O ■ 

• i— i . 



(N. 

>■ 

m ■ 

o, 

^D ■ 
O' 

^t- : 

o, 

o : 

>>: 
^: 

OH- 



X 



Three species collisionless reconnection: Effect of + on magnetotail reconnection 
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The nature of collisionless reconnection in a three-species plasma composed of a heavy species, 
protons, and electrons is examined. Besides the usual two length scales present in two-species 
reconnection, there are two additional larger length scales in the system: one associated with a 
"heavy whistler" which produces a large scale quadrupolar out-of-plane magnetic field, and one 
associated with the "heavy Alfven" wave which can slow the outflow speed and thus the reconnection 
rate. The consequences for reconnection with + present in the magnetotail are discussed. 

PACS numbers: Valid PACS appear here 



Introduction: Recent studies of collisionless recon- 
nection have shown that the disparate masses of the ions 
and electrons lead to a two-scale dissipation region near 
the x-line. The decoupling of the ions from the mag- 
netic field at larger scales than the electrons can lead to 
whistler or kinetic Alfven physics in the dissipation re- 
gion, whose quadratic dispersive characteristics can sub- 
stantially increase the reconnection rate [H 13. [Til fl5l| . A 
Sweet-Parker-like [Hj analysis of this dissipation region 
yields insight into the reconnection rate: Vi n ~ (S/D)cAt, 
where Vi n is the inflow speed, 5 and D are the width and 
length of the dissipation region, and CAt is the Alfven 
speed just upstream from the dissipation region. 

Many plasma systems have heavier species, in addi- 
tion to protons and electrons, which may play an im- 
portant dynamical role: negatively charged dust grains 
in astrophysical molecular clouds and the interstellar 
medium, and + in the Earth's magnetosphere 17]. Pre- 
vious simulations of three-species reconnection have fo- 
cussed on + in the magnetosphere |2l Irl |20| . where the 
number density of + can sometimes exceed that of the 
protons[8|,[l(|. These simulations, however, either did not 
find any effect of + on reconnection due to the small 
system size or did not spatially resolve the reconnection 
boundary layers. 

In this paper we present the first comprehensive study 
of basic three-fluid reconnection showing through theory 
and simulation both the effect of the heavy species on the 
reconnection rate and the hierarchy of scales present in 
the microscale boundary layers. We find that the usual 
two scales associated with collisionless two-fluid recon- 
nection (di — c/ujpi,d e = c/ujp e ) are instead replaced by 
four scales. The inner two scales are associated with a 
light whistler and a light Alfven wave, which are very 
similar to their two-fluid counterparts. At larger scales, 
however, a heavy whistler and heavy Alfven wave occur. 
The heavy whistler can occur on scales much larger than 
a di and thus gives rise to a much wider quadrupolar 



out-of-plane magnetic field signature. Associated with 
this magnetic field are parallel ion Hall currents, the ana- 
logue to the light whistler electron currents. The higher 
+ mass substantially slows the reconnection rate be- 
cause the outflow speed from the x-line is reduced from 
the usual proton Alfven speed, CAi, to the much slower 
heavy Alfven speed, CAt- 

Analytical Analysis: We begin with the general 
three fluid equations, 5j and first ignore electron inertia 
so that E = — V e /c x B — VP e /(n e e). We assume quasi- 
neutrality, n e — m + ZhJih, where n, is the light ion den- 
sity, nh is the heavy species density, and zu is the charge 
number of the heavy species. We ignore the displace- 
ment current, J = (c/47r)V x B. We normalize length 
to di = c^/mi / -\/ Ann^e 2 and time to Q^ 1 — mic/(eBo), 
which gives the following equations: 
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where n e V e = n^V^ + z^n^V^ — J, dV a / dt = (d/dt + 
V a ■ V)V Q , rhu = m h /m l , J = VxB, and P a = T a n a , 
where T a is assumed to be an unchanging spatial con- 
stant (isothermal approximation). 

Although reconnection is a highly nonlinear process, 
much information about its nature can be gleaned from 
a linear analysis, for it is bent field line waves which 
ultimately accelerate the plasma away from the x-line. 
We write each variable as f(x,t) = fo(x) + f e l ( kx -"*). 
Beginning with a uniform Bo with no initial velocities, 
we proceed to linearize Eans. lTl4l and assume that k || Bo 
for simplicity. The sound waves with k || (Vj and V/j) 
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FIG. 1: (below line) The non-ideal length scales present in 
three-fluid reconnection and numbers for typical magnetotail 
lobes with + presents] (n* = 0.05cm -3 , nh/ni — 0.64). 
(above line) Waves and dispersion relations at each scale 
range. 



decouple from the magnetic waves, leaving the following 
dispersion relation: 



± 



Qh z h n h 



ili _ 1,2 ]2 

n t K a s 



Qi K a s 



1 I Qh_ 



Tfc 2 rf^=0, 



(5) 



whe re £lj = eB /(rriic), fl h = ez h B /(m h c), d s = 
diy/rii/n e = c % /ro7/V47rn e e 2 , d h = Cy/rrTh/ ' \J ' \im h z\e 2 , 
n e = rii + Zhnh, and all densities, n, are equilibrium quan- 
tities. This equation is fully general and can apply to any 
plasma with electrons, ions, and a third species. 

Balancing the second and fourth terms and taking the 
limit lo <C ilh S &i and k 2 d\ <C 1 yields the largest 
scale or global Alfven wave: lo = ±kcAt, where CAt — 
B I \J 'Air{mhnh + rmrii). In order for the heavy species to 
slow the global Alfven wave appreciably, it is necessary 
for m h n h > nurn. 

Taking the limit of 51^ 3> £lh, w 3> ^h, and mnTlh 2> 
mirii yields the high frequency dispersion: 
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Zh,n h 



k 2 d: 



k 2 d 2 = 
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For \zh\nh/n e < 1, this equation produces the light 
whistler, the light Alfven and the heavy whistler waves 
shown in Fig. ^ where d\ = & 'rrih / ' (AirnhZ^e 2 ) , and 
CAh = B/y/AimhTrih. The existence of the heavy 
whistler wave has been noted in electron-positron- 
dusty plasmasjlg and electron-proton-dusty plasmas0, 
Il2 |. but was not applied to reconnection. Taking 
k 2 d 2 n e /\zhnh\ 3> 1, and then equating the first and 
second terms yields the light whistler with uo — 
ztk 2 diCAi{ni/n e ), where cai — B / 'v 'Airrmni is the proton 
Alfven speed. Equating the first and third term yields 
the light Alfven wave with w = ±kcAi\/ni/n e . The tran- 
sition between these two waves occurs when k 2 d 2 ~ 1. 
Both of these waves are very similar to their two-species 
analogues except rii has been replaced with n e . 

Taking k 2 d 2 s n e /\z} l n} l \ <§C 1 and equating the first and 
third terms also yields the light Alfven wave. Equat- 
ing the second and third term yields the heavy whistler 
wave with lo = ±k 2 dhCAh- The transition between these 



two waves occurs at k 2 d 2 nin e /{z 2 n\) ~ 1. The heavy 
whistler requires lo ^> Qh so that the heavy species is 
unmoving, but the ion inertia term in Eq. [2]is negligible. 
Thus, the wave is characterized by frozen-in protons and 
electrons that flow together and act as a massless fluid, 
but because rii ^ n e , this net flow is a current. This 
wave transitions to the heavy Alfven wave at kdh ~ 1 
with lo = kcAh in this limit. 

The two scale structure of the dissipation region in 
a collisionless two-fluid plasma[14J (d e and di) has now 
been replaced with four scales: the three scales in Fig. Q] 
plus a very small electron scale S e where the electron 
frozen-in constraint is broken. We did not include S e in 
this calculation to simplify the analysis and because 5 e 
does not appear to substantially modify the reconnection 
rate in well developed Hall mediated reconnection [71 Il3| . 

Simulations: Eans. HMI with Zh = 1 and the same nor- 
malizations were integrated forward in time using F3D, 
a parallel fluid code. The simulation domain is a uni- 
form grid of 2048 x 1024 grid points with the physical 
size L x x L z — 204.8 x 102.4, with periodic boundaries 
at X = ±L x /2 and z — ±L z /2. The initial equilib- 
rium consists of a system size double current sheet with 
B x = £ {tanh[(z + L z /A)/w Q ] - tanh[(z - L z /4)/w ] - 1} 
and wo = 1.5. V/, = initially with n^ = 0.64 every- 
where and Ti = T e = T/, = 0.5. Pressure balance is main- 
tained by setting B 2 /2+{T l +T e )n i = B 2 /2+(T i +T e )n M , 
where Bq = 1.0 and n^ = 1.0 are the values out- 
side the current sheet. A final equilibrium constraint 
is riiViy = J y Ti/(Ti + T e ), which determines Vi y . The 
remainder of the equilibrium current is put into V ey . In 
order to break the frozen-in constraint of the electrons at 
the smallest scales, the term /i4V 4 B has been added to 
the RHS of Eq.EI with ^ 4 = 5 • 10~ 5 . To initialize the 
double tearing mode, x-lines were seeded in both current 
sheets at (x, z) = (±L x /4, =fL z /4) with an initial half is- 
land width w — 0.55. A small amount of random noise 
was added to the initial B and V, of about 10 -3 their 
equilibrium values. 

To examine the effect of the heavy ion mass, We 
ran three simulations with rh = {1, 16, 10 4 } and n/, = 
0.64 in all cases. The first case corresponds to the 
usual two-species reconnection. The second case cor- 
responds to reconnection in the presence of + , with 
{d s ,d iy /nin e /(z h n h ) 2 ,d h } = {0.8,2,5}. In the third 
case, the 3 length scales are {0.8,2,125}, so that the 
heavy ions form an immovable background and the global 
scales are controlled by the heavy whistler. 

The reconnection rates of the three simulations versus 
time are shown in Fig.|2Jt. The rhh = 1 case clearly shows 
a substantially larger reconnection rate. The two cases 
with heavy ions show very similar reconnection rates, but 
the largest rhh shows a large decrease in its reconnection 
rate around t = 600, while the rhh = 16 case keeps a 
steady rate. The heavy whistler velocity has a k de- 
pendence, V ~ kdhCAh- As reconnection proceeds in a 
system and the island width w gets larger and larger, 
the effective k ~ 1/w for the reconnection process de- 
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FIG. 2: (Color online) (a) Reconnection rates, (b)-(d) B y 
with proton flows, (b) mj, = 1 and t = 500, (c) mj, = 16 and 
t = 650, (d) rh h = 10 4 and t = 650. 



creases. Because the heavy whistler is mediating global 
convection in the rhh — 10 4 case, as the global convec- 
tion scale length increases, the global convection velocity 
must decrease, throttling the reconnection rate. 

The reconnection generates very different signatures 
for the different rhh ■ Figs. |2>d show the out-of-plane B y 
generated from the reconnection and the proton flow vec- 
tors. The x-line is located close to (x, z) = (10,-25.6) 
in all three cases. Only a small fraction of the total sim- 
ulation is shown. The rhh = 1 case shows the usual 
quadrupolar structure generated by frozen-in electron 
flow[9|. For x > 45, the clean quadrupolar signature 
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case, a cut along z at x = 20.0 : (a) B y 
The vertical dotted line is the symmetry 



begins to change to a more complicated structure with 
both positive and negative B y due to the finite system 
size. Because Vi X is maximum at about x — 45, and for 
greater x the slow-down of V% x causes a compression of 
B z , and the resulting J y generates a B y signature of the 
opposite sign. 

The rhh = 16 case (Fig. Eh) shows signatures of both 
the light and heavy whistler. For x < 40 a narrow band of 
positive By associated with the light whistler is present. 
A cut of B y , Vi X , and Vhx at x = 20 is shown in Fig. 
This B y spike has a main length scale of about a di, 
which is roughly consistent with the light whistler cut- 
off scale of d s = 0.8 for this simulation. Like the usual 
two-species whistler, the in plane current generating this 
By is due to counterstreaming parallel electron beams 
upstream and downstream of the B y perturbation. There 
is a long tail of B y upstream of the spike (z > —23.5) in 
Fig. Et, though, which is not present in the two-species 
case. The proton outflow shows a peak on the symmetry 
axis like the two-fluid case, and its velocity is much larger 
than the + velocity. 

The quadrupolar B y becomes dominated by the heavy 
whistler for x > 40 in the m^ = 16 case. The B y signa- 
ture broadens out substantially because dh — 5 for this 
simulation, and the current which generates it is carried 
by both the ions and electrons. Fig. 0Ji shows compar- 
ison slices for the rhh = 1 and 16 cases at x = 55.0. 
The main positive B y spike is substantially wider in 
the rhh = 16 case, although it is not 5 times wider as 
might be expected from a comparison of dh to di. The 
x-velocities reveal another key signature, as shown in 
Fig. 0Jd for rhh — 16. The parallel ion flows from the 
heavy whistler associated with B y lead to a negative Vi X 
at about z = —17. Also, the ion flow no longer is maxi- 
mum at the symmetry line, but instead peaks off axis at 
around z = —20.5. On the symmetry line, Vhx is some- 
what larger than V, x - Vi x is about 4 times slower in the 
rhh = 16 case than in the rhh = 1 case. The off axis 
peak of V,x and the substantial negative Vi X (about 1/3 
of maximum ion outflow) do not occur unless the heavy 
whistler is active. 

In the case with rhh = 10 4 , the heavy ions are immov- 
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FIG. 4: (a) Slice of B y along z at x = 55, (b) slice of x- 
velocities along z at x = 55 for rhh = 16, (c) slice of x- 
velocities along x at z = —25.6 for rhh = 16. 



able and the heavy whistler is dominant at the global 
scales in the simulation as seen in Fig. EJi. The main 
peak of By is quite wide, and there is a nonzero B y out 
to global scales. The parallel ion flows which flow with 
the Hall electron currents are much stronger in this case. 
The multiscale structure of the dissipation region is 
demonstrated for the + case in a cut of the outflows 
away from the x-line, as shown in Fig. 0J:, which is a 
cut along x through the center of the current sheet. For 
x < 21, the behavior of the flows is v ery similar to that 
seen in two-species reconnectionflol fl6) . In the light 
whistler region, the electrons accelerate to speeds much 
faster than the ions and exceeding the relevant Alfven 
speed. The electrons cannot maintain this high velocity 
and are forced to decelerate to a speed comparable to the 
protons when they reach the light Alfven region. Inside 
of this Alfven region, the protons reach their maximum 
velocity. The protons, in an analogy to the electrons 
in the light whistler region, cannot flow this speed in- 
definitely, and are forced to slow down inside the heavy 



whistler region. Finally, their speed becomes compara- 
ble to the + outflow in the heavy Alfven region. Inside 
the heavy whistler region at x = 29, however, Vi x drops 
nearly to zero, below Vh x - This behavior also occurs in 
the rhh = 10 4 case at about x = 18, as seen in Fig. [5Ji. 
The + outflows behave like the proton outflows in the 
two- fluid case, gradually accelerating and finally reaching 
their maximum velocity in the outermost Alfven region. 

Discussion As discussed in the introduction, a Sweet- 
Parker like analysis of the dissipation region yields 
Vin ~ {5/D)cAt- For the rh h = {1,16} cases, we 
would expect the outflow speeds to differ by a factor 
of y/(minn6 + ThhTih) / (niirni) — 2.6 between the two 
cases, where nn = 1.64, n,i6 = 1.0, and rih = 0.64. The 
maximum outflow in the two-species case is about 1.0, 
while for the case with + is .35, giving a ratio of 2.9, 
quite close to what is expected. If S/D stayed the same 
between the two simulations, the + case would be ex- 
pected to reconnect nearly 3 times slower than the two- 
species case. In Fig. |2K, however, the two cases asymp- 
tote to approximately steady-state rates that differ by 
around 1.5, leaving a factor of about 2 unaccounted for. 
A rough estimation of the scaling of S/D between the sim- 
ulations may be possible by examining the angle, 0, that 
the By signature makes with the z — —25.6 symmetry 
line. Presumably S/D ~ tan6>. This gives tan#; = 0.08 
and tan$o+ — -15, which sheds light on the factor of two 
difference. A more careful determination of S/D as well 
as a scaling study with very large system sizes will be 
necessary to determine if this change in S/D is robust. 

These results imply that a substorm occurring with a 
high enough + density (mhUh 3> W,rij) will have slower 
outflows and a reduced reconnection rate normalized to 
the upstream proton Alfven speed. With all else being 
equal, this implies that the expansion phase of substorms 
will take longer to occur or will reconnect less lobe field 
in the same amount of time. However, substantial + 
populations tend to occur during times of increased ge- 
omagnetic activity. The magnetotail equilibrium, and 
thus the lobe magnetic fields and density, may be modi- 
fied substantially during these periods, which may offset 
or even overpower the reduction in the reconnection rate. 
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